Preprint typeset in JHEP style - HYPER VERSION 



April 30, 2010, KEK-CP-234, RBRC-842 



The epsilon expansion at next-to-next-to-leading 
order with small imaginary chemical potential 



Christoph Lehner"^, Shoji Hashimoto'^, and Tilo Wettig" 

Institute for Theoretical Physics, University of Regensburg, 93040 Regensburg, Germany 
^RIKEN/BNL Research Center, Brookhav en National Laboratory, Upton, NY- 1197 3, USA 
'^High Energy Accelerator Research Organization (KEK), Tsukuba 305-080, Japan 

Email: clehner@ quark . phy . bnl . gov, shoji . hashimotoSkek . jp, 

t ilo . wett igSphysik . uni-regensburg . de 

Abstract: We discuss chiral perturbation theory for two and three quark flavors in the epsilon 
expansion at next-to-next-to-leading order (NNLO) including a small imaginary chemical potential. 
We calculate finite- volume corrections to the low-energy constants S and F and determine the non- 
universal modifications of the theory, i.e., modifications that cannot be mapped to random matrix 
theory (RMT). In the special case of two quark flavors in an asymmetric box we discuss how to 
minimize the finite- volume corrections and non-universal modifications by an optimal choice of the 
lattice geometry. Furthermore we provide a detailed calculation of a special version of the massless 
sunset diagram at finite volume. 
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1. Introduction 

At low energies the theory of quantum chromodynamics (QCD) can be described by a chiral effec- 
tive theory. If the theory is confined to a finite volume and considered for small quark masses, the 
e-regime power counting applies and replaces the standard p-regime power counting, see Ref. [1]. 
The con^esponding systematic expansion is called e-expansion. 

At leading order (LO) in the e-expansion the theory becomes zero-dimensional and is there- 
fore described by chiral RMT [2], see [3, 4] for reviews. The dimensionless quantities of RMT 
ai^e mapped to the dimensionful quantities of the chiral effective theory using the LO low-energy 
constants (LECs) S and F, see, e.g., Ref. [5]. In this way S and F, which are of great phenomeno- 
logical importance, can be obtained from lattice QCD simulations in the e-regime by a fit to RMT 
predictions. While S can be determined rather easily, e.g., from the distribution of the small Dirac 
eigenvalues, the extraction of F is more complicated and can be done, e.g., by the inclusion of a 
suitable chemical potential [6, 7] or by using twisted boundary conditions [8] . 

Lattice QCD simulations in the e-regime using an exactly chiral lattice fermion formulation 
are feasible, and recently such a set of lattice configurations, which also includes the effect of 
dynamical quarks, was generated by the JLQCD and TWQCD collaborations [9, 10, 11]. In an 
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upcoming publication [12] we further investigate the eigenvalue spectrum of the Dirac operator on 
these configurations. 

For typical volumes of state-of-the-art lattice QCD simulations, finite-volume corrections to 
the RMT predictions cannot be neglected. These corrections can be calculated in higher orders of 
the e-expansion. At next-to-leading order (NLO) the analytical results of RMT still apply, but the 
mapping of RMT to chiral perturbation theory is modified, i.e., the LECs S and F are replaced by 
finite-volume effective LECs Sgff and Feff, see Refs. [13, 14, 15]. At NNLO, chiral perturbation 
theory can no longer be mapped to RMT, i.e., non-universal modifications appear. These non- 
universal modifications determine the systematic errors in fits of lattice data to RMT predictions. 

In this paper we calculate the finite-volume corrections at NNLO in the e-expansion in Eu- 
clidean space-time. We allow for nonzero imaginary chemical potential and consider its contribu- 
tion to leading order. In this way we extend results of Hansen [16], which were obtained without 
chemical potential. The paper is structured as follows. In Sec. 2 we calculate the finite-volume ef- 
fective action and the corresponding finite-volume effective low-energy and high-energy constants 
(HECs) at NNLO. In Sec. 3 we discuss how to minimize the systematic deviations from RMT as 
well as the finite-volume coiTcctions to S and F in the special case of two quark flavors on an 
asymmetric lattice by an optimal choice of the lattice geometry. We conclude in Sec. 4 and provide 
a detailed calculation of a special version of the massless sunset diagram at finite volume in App. A. 

2. The finite-volume effective theory at NNLO 

In this section we discuss the finite-volume effective theory at NNLO in the e-expansion. To this 
end we give the bare Lagrangian of the p-expansion at NLO and the LO contribution of the invariant 
integral measure in Sec. 2.1. In Sec. 2.2 we define the finite- volume effective action, and in Sec. 2.3 
we calculate the finite- volume effective LECs and HECs. The LECs Li, . . . , Lg and the HEC H2 
that appear in the bare Lagrangian of the p-expansion at NLO are scale-dependent. We renormalize 
the theory in Sec. 2.3 and confirm that the scale dependence of Li, . . . , Lg and H2 is the same as 
in the ordinary p-expansion [17]. 

2.1 The bare Lagrangian and the path-integral measure 

We parametrize the Nambu-Goldstone manifold of chiral perturbation theory with Nf quark flavors 



where ^ is a complex matrix in flavor space of dimension Nj with ^ = and Tr ^ = 0. The 
constant mode is separated in Uq, and thus J d'^x£^{x) = 0, see Ref. [1]. For nonzero imaginary 
chemical potential the LO Lagrangian of the effective theory in Euclidean space-time is given by 



by 



U {x) = Uq exp 



^^/2 



(2.1) 



F 




(2.2) 



with 



VpU{x) = dpU{x) - i5po[C, U{x)] 



(2.3) 
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where M = diag(mi, . . . ,mj^^.), ruf is the mass of quark flavor f,C = diag(/ii, . . . ,iJ.Nf), 
and i^f is the imaginary chemical potential of quark flavor /, see, e.g., Ref. [14]. We consider 
very small quark masses such that the Compton wavelength of the pions, given by the Gell-Mann- 
Oakes-Renner relation 

ml = , (2.4) 

exceeds the size of the space-time box. This defines the e-regime of QCD that was first discussed 
in Ref. [1]. The corresponding e-regime power counting [1] is defined by 

y ~ e"^ , M ~ e^ C ~ e^ 5p ~ e , ^(x) ~ e (2.5) 

which gives a consistent perturbative expansion in ~ l/F^-y/F, see, e.g., Ref. [15]. The e- 
expansion is applicable if the quantities 

m/FS and njP'^V (2.6) 

are not larger than 0(1) and the volume is sufficiently large, i.e., 

AttF'^VV > 1 , (2.7) 

see also the quantitative discussion in Sec. 3. 

In order to include all NNLO contributions in the e-expansion we need to include the NLO 
Lagrangian of the p-expansion which for Nf = 3 is given by 

£2 = -LiiTiiV^U^^Vf'U)]'^ - L2Tr{V^U^^VuU)TiiV^U~^V''U) 
j L5 Tr[V^C/-^V^?7(M[/-i + M^U)] 

j Ls Tr{MU'^MU~^ + M^UM^U) -[^) ^2 Tr(MtM) , (2.8) 

where H2 is a high-energy constant (HEC) corresponding to a contact term that is needed in the 
renormalization of one-loop graphs, see Refs. [16] and [17]. The field-strength tensors are not 
included since they vanish in the case of a constant vector source [17], of which an imaginary 
chemical potential is a special case, and therefore the LECs Lg and Liq and the HEC Hi do not 
appeai\ Note that in the case of Nj = 2 not all of the terms in £2 are independent and the LECs 
Li , . . . , Lio can be mapped to the LECs Zi , . . . , ^7 as described in Ref. [18]. In the case of Nf > 3 
one needs to include additional terms in Eq. (2.8), see, e.g., Ref. [18]. 

The invariant measure relevant to NNLO in the e-expansion is given by 

d[U] = d[UoW] (1 - ^ I d^xTiUxf]^ , (2.9) 

see Ref. [16], where d[Uo] is the invariant measure of the constant mode and d[^] is the flat measure 
of the fields ^. 
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2.2 The finite-volume effective action 

The volume dependence of the theory is contained entirely in the propagators of the fields ^ [19], 
and therefore we can obtain a finite-volume effective action in terms of the constant mode Uq 
by averaging over the fluctuations in ^. At NLO this leads to a finite-volume effective action 
that differs from the LO (RMT) result by finite-volume connections to the LECs S and F, see 
Refs. [13, 14, 15]. We now perform the expansion of the partition function 




(2.10) 



in terms of fields ^ to NNLO and average over the fields using computer algebra, resulting in 

Z = j d[C/o]e-^=«. (2.11) 

The finite-volume effective action S'eff contains invariant terms including the constant mode Uq, the 
mass matrix M, and the chemical potential matrix C. At NNLO in the e-expansion and to leading 
order in the imaginary chemical potential, S'eff can be written as 

S'eff = Tr(Mt[/o + U^'M) - ^ Tr{CU^'CUo) 

+ Ti^iVFf Tr(C)[Tr(C/o{Mt, C}) + Tt{Uq^{C, M})] 
+ T2^{VFf Tr({Mt, C}UoC + {C, M}CUq^ 

+ {Uo, C}U^^CUoM^ + CUo{C, U^'}MUo^) 
+ Ts^iVFf Tr{UQ^CUoC + C^) Tt{MU^^ + M^Uq) 
+ T^J:{VFf Ti{Uq^CUoC - C^) Tt{MU^^ + A'PUo) 
+ T5^{VF)^ Tr([Mt, C]UoC + [C, M]CUq^ 

+ [Uo, CpQ^CUoM^ + CUo[C, U^^]MUq^) 
+ T6(F$])2[Tr(MC/o-i + Mt[/o)]2 
+ T7(yS)2[Tr(MC/o-i - M^Uo)f 
+ T^{VT.f[TT{MU^^MU^^) + Tv{M^UoM^Uq)] 

+ 'HiVF^TY{C'^) + n2{VT.fTY{M^M)+'H^VF'^{TTCf, (2.12) 

where {^4, B} = AB + BA and Seff, Fgff, and Tf, . . . , Tg are finite- volume effective LECs that 
will be given in Sec. 2.3. Note that the terms corresponding to T-Li, 'H2, and T-L^, do not couple to 
Uq, and therefore the T-Li can be viewed as finite-volume effective HECs. The terms proportional to 
the T-Li only contain external sources and are therefore not relevant for low-energy phenomenology, 
but they are needed for the computation of operator expectation values, see Ref. [20]. They are also 
relevant to the renormalization of the coupling constants discussed in Sec. 2.3. The finite-volume 
effective HECs T-Li are also given in Sec. 2.3. 

Unlike the LECs Li, . . . , Lg and the HEC H2, the finite-volume effective LECs Tf, . . . , Tg 
and HECs T-Li, 1-12, and T-Ls are finite and depend on the volume. Specifically, we show in Sec. 2.3 

'We use a C++ library for tensor algebra developed by one of the authors. 
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that 

Seff,Feff,^i = 0(e°), n^ = 0{e^), n2,Ti,... ,Ts = 0{e^), (2.13) 

i.e., 712, Tis, and Ti, . . . , Ts vanish in the infinite-volume limit. 

The terms corresponding to the Tj cannot be mapped to RMT, see, e.g., Ref. [5]. Therefore 
their magnitude quantifies the systematic deviations from RMT at finite volume. We will return to 
this point at the end of Sec. 3. 

2.3 The finite-volume effective low-energy and high-energy constants 

In the following we express the finite-volume effective LECs Seff, F^g, and Ti, . . ., Tg and the 
finite- volume effective HECs T-Li, 'H2, and T-L^ in terms of shape coefficients Pi, . . . ,Pq defined be- 
low. The resulting expressions for the effective LECs and HECs are given in terms of the massless 
finite-volume propagator in dimensional regularization, 

1 ikx 

where the sum is over all nonzero momenta, see Ref. [19]. We use the identity 

dlG{x)l^, = ^ (2.15) 
and finally express the result in terms of the shape coefficients 

Pi = VdlG{{)) , P2 = VVG{0) , 

p^ = ^/V j d'^x[diG{x)]G{x) , j '^'^^^(^)^ ' 

P5 = y d'^x <fy[dlG{x + y)]G{x)G{y) , P^ = ^ J d''x[diG{x)]G{xf , (2.16) 

where d is the number of space-time dimensions. We use conservation of momentum, by which all 
one-loop propagators can be related to [19] 

where r G R and r(r) is the Gamma function [21]. The partial derivatives can be evaluated using 

where Ly is the length of the space-time box in dimension v = 0, 1, 2, 3 and denotes the 
partial derivative w.r.t. Ly. Note that in Eq. (2.18) no sum over v is implied. The shape coefficients 
Pi , . . . , P5 contain only one-loop propagators and are given by 

Pi = ~di^fioGo), P2 = ^Gu 

2" 

P^ = -\dj^^{UG2). (2.19) 
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For convenience we state the result of Ref. [19] explicitly as 

1 



Gr = lim 



m— >0 



(4vr)'^/2 



T{r - d/2){m^)'^l'^~'' + gr 



r(r) 



Vgo = /3o + Pim'^W + ^/SsmV - log(mW) 



+ 



9r+l 



Vrn^ 
2(4^ 

dgr 
d{m?) 



log(m vV) 



(2.20) 



with shape coefficients /3„ given in Eq. (B.14) of Ref. [19]. We express Go, Gi, and G2 in terms 
of /3o, /?!, and /32 and find 



P5 = ->4-^Lo5^/2 



-/3i 



-2A + /32 + 



logV^ 
(47r)2 



4- ' 4 -1.0'- (i6vr)2 ' 
where we borrow the definition of A from Ref. [17], 



(2.21) 



A 



d-4 



1 



1 



(47r)2 \ d-A 2 



(47r)2 \ d-A 2 



[1 + r'(i) + iog(4^ 



[1 + r'(l) + log(/i2) + log(47r)] 



-)]} 



(2.22) 



We explicitly include the dependence on the scale ^, which we define with mass dimension one. 
Note that the logarithms of dimensionful quantities in Eqs. (2.21) and (2.22) can always be com- 
bined to logarithms of dimensionless quantities. The shape coefficient Pq is calculated in App. A, 
see Eq. (A.41). It contains a special version of the massless sunset diagram at finite volume. 

In the following we state the resulting expressions for Sgff, Fgff, Ti, . . . , Tg, T-Li, 7^2, and T-L-j,. 
The finite- volume effective chii^al condensate is given by 



■ 1 

+ 



F^VV 
Pa 



2'' "I ' F^V 

which agrees with Eqs. (22) and (23) of Ref. [16]. The finite-volume effective LEC 



(2.23) 



F2 ^ F^VV 



2Nf + 2N'j^^ + 2N^^^ + N'j 



fWv ■ --fFW^^'^fm^^'^fFW 



F^V 



+ ^ [2Li + N]L2 + {Nf - 2NJ^)L, 



F^V 



(2.24) 
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and the finite-volume effective HEC 



+ ^ [(iV| - l)Li + L2 + (iV/ - Nj')Ls] 

+ 1^ [2Li + N]L2 + {Nf - 2Ny')L3] (2.25) 

contain the contribution of the two-loop propagator in Pq. The other finite- volume effective LECs 
and HECs are given by 

1 P4 + 1 AT -Y^ 

^1 = 2 W ' T2 = --iV/Ti, 

^ V „ 1 P4 4L4 



2 4 F^F ' 

1 2iv5 1 , 2\ 4Zyg 



and 



(P2 + 2P3) iVf (Pn^ + 4P2) (P4 + 4P5) P2P3 

Note that the shape coefficients P4, P5, and Pg as well as H2 and the Lj are divergent and need 
to be renormalized. We separate their scale dependence as 

P4 = PI-2A, P, = P^ + ^X, p^ = pr + lx-^P,X, 

L, = L\ + TiX with i = 1, . . . , 8 , H2 = H'2+A2X, (2.28) 

where the quantities with superscript r are finite. For Nf = 3 the divergences in Eqs. (2.23) - 
(2.27) can be absorbed if and only if we choose 

1 3 11 

^^-8' ^^-8' ^^-144' 

rr = 0, rs = ^, A2 = ^, (2.29) 

and 

30 16 T 

— = iQTi + 2r2 + — = 2ri + 9r2 + -rg . (2.30) 

lb 3 3 

The coefficients r4, . . . , Fg and A2 are equal to the coefficients obtained in the one-loop expansion 
in the p power counting, see Ref. [17]. The renormalization conditions of Eq. (2.30) for Fi, F2, F3 
are also compatible with the result of Ref. [17], 

Fi = — , F2 = — , F3 = 0. (2.31) 
^ 32 ' 16 ' 
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Note that the divergences in H.^, Ti, T2, and T3 cancel independently of the choice of Fj. 

To summarize, we obtain finite expressions for the finite-volume effective LECs and HECs in 
Eqs. (2.23) - (2.27) if we replace P4, P5, Pq, H2 and the Li by their corresponding renormalized 
parts with superscript r. Note that the dependence on the scale ^ drops out of the final results 
for the finite-volume effective LECs and HECs. The renormalization in the case of Nf = 2 is 
discussed in Sec. 3. 



3. Optimal geometries: Two quark flavors in an asymmetric box 

In the following we discuss the finite-volume coiTcctions to S and F and the coefficients of the 
non-universal terms at NNLO in the e-expansion. We explicitly consider the case of Nj = 2 and 
an asymmetric box with lattice geometries 

{a-x) Lo = xL, Li = L2 = L3 = L, 

{bx) L3 = xL, Lo = Li = L2 = L, (3.1) 

where x G {1, 3/2, 2, 3, 4}. The three-flavor coupling constants Li, . . . , L5 can be related to the 
two-flavor coupling constants li, I2, and I4 by 

h = 4Li + 2L3 , I2 = 4L2 , U = 8L4 + 4L5 , (3.2) 

see, e.g., Eqs. (3.15) and (3.16) of Ref. [18]. Therefore 

Seff ^ 3P2 _ 3P| 3^ 3Z^ 

S 2F'^\/V 8F^V F^V F^V ' 



iffF^ 4P2 4P3 8P2P3 + 8PI+4PI 8P4 + I6P5+4P6 



and 



12^1 + 4/2 ^1(8/1 + 16/2) 
In Ref. [22] the scale dependence of the coupling constants li with i = 1, . . . , 7 is separated as 

/^ = /[+7^A, (3.5) 

where 

1 2 

71 = 3 , ^2 = 3 > 74 = 2 . (3.6) 

It is straightforward to check that the divergences in Eqs. (3.3) and (3.4) cancel with this set of 71, 
72, and 74. We therefore obtain finite results for Egff and Feff if we replace ^4,^5, Pe and /i, ^2, ^4 
in Eqs. (3.3) and (3.4) by their corresponding renormalized parts with superscript r. 

The shape coefficients Pi, P2-,Pz and the renormalized shape coefficients PJ , , Pg at scale 
/I = are given in Table 1 for geometries (clx) and ^x)- The details of the calculation of Pg 

are given in App. A. 
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P\ 


P2 


-P3 


PI 


pr 


Pe 


(ai) 


0.249999 


—0.140461 


0.035115 


-0.020305 


0.004285 


-0.0106(1) 


(03/2) 


-0.251470 


-0.123339 


—0.014414 


-0.018457 


0.001121 


0.00940(4) 


(02) 


-0.674807 


-0.083601 


-0.076021 


-0.012954 


-0.005093 


0.02487(2) 


(as) 


1 c; 1 o^^ 1 n 
— i.OlzDiU 




—U./oi^l 1 


n ni 01 p^n 


— U.Uoiibo 


n HQi af 1 ^ 
U.Uoiy^i ) 


(04) 


-2.350148 


0.215097 


-0.440882 


0.062628 


-0.082324 


-0.0214(3) 


(^3/2) 


0.417156 


-0.123339 


0.045917 


-0.018457 


0.004723 


-0.01513(3) 


(^2) 


0.558268 


-0.083601 


0.053207 


-0.012954 


0.004960 


-0.01577(7) 




0.837535 


0.041942 


0.065178 


0.012150 


0.005282 


-0.0080(1) 


(^4) 


1.116713 


0.215097 


0.075261 


0.062628 


0.005510 


0.0117(1) 



Table 1: Shape coefficients Pi, P2, ^3 and renormalized shape coefficients P^, P^, Pq at scale ji ^ V 
for geometries (fla,) and (bx). The error in the last column is due to the extrapolation described in App. A. 3. 







(^3/2) 


(«2) 


(«3) 


(04) 




1.3455 


1.2477 


1.1454 


0.9404 


0.7355 




1.39(1) 


1.288(7) 


1.202(5) 


1.047(3) 


0.906(3) 




1.3004 


1.3182 


1.3192 


1.3193 


1.3193 


^NNLO/^ 


1.279(9) 


1.305(4) 


1.306(2) 


1.292(1) 


1.261(2) 






(^3/2) 




(bs) 


(^4) 


pNLO/p 




1.1894 


1.06816 


0.7710 


0.2186 


pN^LO/p 




1.182(8) 


1.092(7) 


0.959(6) 


0.919(5) 



Table 2: Finite-volume corrections to E and F at NLO and at NNLO for geometries (a^) and (bx) at 
Tn^VV = 1, F = 90 MeV, and L = 1.71 fm. The error in is due to the uncertainty in [4, the error in 

F^^'^ is due to the uncertainty in Zi, ^2, and Pq. 

The renormalized coupling constants can be related to scale-independent constants li by 

i: = ^^[li + logiml/t,^)], (3.7) 

where m.,^ is the mass of the pion, see Ref. [22], and 

[1 = -0.4 ±0.6, [2 = 4.3 ±0.1, [4 = 4.4 ±0.2, (3.8) 
see Ref. [20]. Therefore 

i: = ^^[li + logimlVv)] (3.9) 

at scale 1^1 = Note again that the finite-volume corrections to S and F are independent of 

the choice of scale fi. 

In Table 2 we give explicit values for the finite-volume corrections to S and F at NLO and 
NNLO for geometries {a^) and (bx) with m\\/V = 1, F = 90 MeV, and L = 1.71 fm, which 
roughly corresponds to the values of the JLQCD lattice simulations. Note that the e-expansion 
converges well for this set of parameters as long as the asymmetry of the lattice is not too strong 
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Figure 1: Finite-volume corrections to S and F at NNLO for geometries (a^) and (b^r) at m'^^/V = 1, 
F = 90 MeV, and L = 1.71 fm. 



(convergence is worst in geometry (64)). We calculate the NLO result by discarding terms of order 
in Eqs. (3.3) and (3.4). Note that for the same value of x, Seff is independent of the choice 
of lattice geometry (ax) or (bx). In Figure 1 we visualize the NNLO results of Table 2. We confirm 
the picture obtained in Ref. [15] at NLO that the finite- volume connections to F can be largely 
reduced by an asymmetric lattice geometry with one appropriately large spatial dimension instead 
of one large temporal dimension. 

We now turn to the non-universal terms that cannot be mapped to RMT. It follows from 
Eqs. (2.21) and (2.26) that the coefficients T4, . . . , Tg are independent of the choice of lattice 
geometry (a^) or (bx) for the same value of x. The coefficients Ti, T2, and T3, however, are 
affected by the choice of lattice geometry (qx) or (bx), and we have 

Tl,T2,T3 0cPI + 4P5^ (3.10) 

We give values for PI + AP^ in Figure 2 at scale = V~^^^ for lattice geometries (ux) and (bx)- 
Note that the non-universal contribution of Ti, T2, and T3 is reduced significantly in lattice geom- 
etry {bx) compared to lattice geometry (a^,) for the same value of x. In an upcoming publication 
[12] we perform the corresponding lattice simulation for x = 2 and show numerically that the 
systematic deviations from RMT are indeed smaller for lattice geometry (62) compared to lattice 
geometry (02). 
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• Geometries {ax) 

♦ Geometries {bx) 
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Figure 2: Linear combination of shape coefficients + iP^ for different geometries at scale 11 ^ V ^1^. 
We divide by x since V = xL^ in Eq. (2.26). 
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4. Conclusions 

We discussed the e-expansion at NNLO and determined the finite-volume effective action and the 
finite-volume effective LECs and HECs to this order. In the special case of two dynamical quarks 
confined to an asymmetric box we have confirmed the picture obtained at NLO that finite-volume 
corrections to the LECs S and F can be significantly reduced by choosing one appropriately large 
spatial dimension instead of a large temporal dimension, see Ref. [15]. Furthermore, we have 
shown that the systematic deviations from random matrix theory can also be reduced in the setup 
with one large spatial dimension. This implies that in order to determine the LECs S and F from 
eigenvalue correlation functions, as suggested in Refs. [6, 7] and performed in a pilot study in 
Ref. [23], one should choose an asymmetric lattice with one large spatial dimension. This will be 
demonstrated numerically in an upcoming publication [12]. 

We would like to add that even though we did not explicitly perform our calculations in a 
partially quenched setup, it is straightforward to extend our results to the partially quenched case, 
see, e.g., Ref. [15] for a discussion at NLO, and we expect our findings to be unmodified by the 
presence of valence quarks. 
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A. The massless sunset diagram at finite volume 

In this section we calculate the two-loop contribution defined by 



where the sum is over all nonzero momenta k and p with p + k ^ 0, and ko is the temporal 
component of the momentum vector k. We first express the propagators without constant mode as 
the limit of ordinary, massive propagators. 




(A.1) 



lim Pe{m^) = lim \ {m'^) + {m^ 



)] 



(A.2) 



with 




(A.3) 



The terms P^{m?) and Pg^(m^) are calculated separately in the following. 
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A.l The term Pg°(m2) 

We partition the term PQ{m'^) into its infinite-volume part and the finite-volume propagator gi 
defined in Eq. (2.20). We find 



E 



where 



gi{w?) 



1 /3i ?n2log(m2/V^) 



Vm? 



V 



(4vr) 



771^/32 + C'(m^). 



(A.4) 



(A.5) 



Therefore we can express P^{m?) in terms of jSi, (^2, and A as 
1 



Pe^K) = 2A 



2m2\/y 



A.2 The term Pg^(m2) 
The second term 



1 



V ^ k'^ + [p^ + m^) [(p + /c)2 + m^] 



is more involved. We use Poisson's summation formula 



n 



(A.6) 



(A.7) 



(A. 8) 



and write 



Pl{m') = -VY. j 



(2^(2^ 



exp 



\^Lj{rjkj + SjPj) 



(p^ + m?) [{p + k)"^ + m?] {k"^ + m?) 



where the sum is over r, s G Z"^. We partition the sum over r and s into 

[A) r = OAs = 0, 

{B) r/OAs = 0, 

(C) r = OAs/0, 

(£>) r/OAs/OAs = r, 

(^) r/OAs/OAs/r. 



(A.9) 



(A. 10) 



Part (A) is given by the infinite- volume sunset diagram, see Ref. [24], which scales with Vm'^ and 
therefore vanishes in the massless limit. The parts PQ^{m?), . . . , P^^{m?) are calculated in the 
following. 
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A.2.1 The term Pe^^(m2) 

Along the Unes of Eqs. (A. 10) and (A. 11) of Ref. [25] we separate 



(A. 11) 



with 



prim') 



p 



1B2, 



m 



d'^k d'^p kl 1 



{2ttY {27tY + m 
1 



exp I 



The term Pg^^^(m2) contains the ultraviolet divergence and can be calculated explicitly, 



(A. 12) 



-V 



d'^p 



{2TT)d (p2 + ^2)2 



E 



d'^k 



(A. 13) 



(47r)2 

where Pi is the one-loop shape coefficient defined in Eq. (2.16). The term Pq^^ (m^ ) is finite. After 

6 



a tedious but straightforward calculation along the lines of Ref. [19] we can express P^^'^{m?) as 



,1_B2 



1 ■ \ 

— / dxdydz K{x,y,z) exp -{x + y + z)- 



47r 



(A. 14) 



with 



K{x,y,z) 



1 



2(x + y) 



(L2/Vy)[2ro(x + y)]27rn 



{yz + xy + xz) 



X exp 



E 



1 



(x + ?/)2z3 



{Ly^V)Tj{x + y)n 
{yz + xy + xz) 

2 - {Ll/Vv)^7:rl/z 



exp 



„2n 



(A. 15) 



This expression is suitable for a numerical evaluation of PQ^'^{m?) if we perform the integral over 
x, y and z in spherical coordinates. In Sec. A.2.4 we discuss how to efficiently calculate infinite 
sums such as the sums over ro, . . . , rs with r 7^ in Eq. (A. 14). 

A.2.2 The term Pgi'^(m2) 

The method used to separate the divergent part of Pg^^(?n2) does not work for the integral over k 
since it has a power divergence. Nevertheless, we can calculate the divergent sub-diagram 



Ifj,u{m,p) 



d'^k 



kfxki/ 



(27r)'^ (fc2 + 7v?)[{p + A;)2 + m^] 



(A. 16) 
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explicitly. The result is given by [26] 



^ , , , fm^ + x(l — x)p'^] log[m^ + x(l - x)p'^] 
I,iu{m,p) = g^u I ax ofA^\2 



2(47r)^ 

PfiPu J dxj^^ [1 + log[m^ + x(l - x)p'^]] 



9fiiy^[ ^P^ +mA - '^Xp/.Pu ■ (A. 17) 



We can thus separate the divergent part of 



Pg (jn) = -V^j ^^2^^exp(^i^L,s,p,J , (A.18) 



which is given by 



)]uv = -^E y '-^w[^Y.L,s,p, 



5, 2_ 2 



A--AP1 + OK). (A. 19) 

b a 

In the calculation of [Pg^'^(m^)]uv we used the identities 

E / exp(^z^L,s,pjj =0, 

E / (2!).p2f^2 E ^-..P.) = -Pi + O(m^) . (A.20) 

Note that the first two identities hold for arbitrary d. Thus there is no finite contribution from the 
product of these integrals with A. 

The finite contributions to PQ^{m?) are given by 



where 
with 



Fi{m^ ,p^) = [m^ + x(l — x)p^] log[m^ + x{l — x)p^] , 

J"2(m^/) = 2x^(1 + log[m2 + -x)/]) . (A.23) 
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<' i\/v\ + TP? Ix(\ — x) 



Figure 3: The complex plane of po- 



We define L\ = LiSi and rotate the coordinate system of p such that 



J'Jrr?) 



dx 



■ exp(iL'*po) 



with = [Z]n=o(^n)^]^''^- After differentiating w.r.t. so we find 



(2vr)^ 



with 



and 



2 2,2 



(A.24) 



(A.25) 



(A.26) 



(A.27) 

In Figure 3 we sketch the structure of the integrand of Eq. (A.25) in the complex plane. There 



are two poles at pQ = ±i^p'j_ + m? and a branch cut due to the logarithms in Ti{m? ,p^) and 
F2{'m? tP^)- We can close the integration contour in the upper half-plane and find 



j^K) = [j^K)]p + [j^K)] 



(A.28) 



where [J's{m'^)\p is the contribution of the pole and [J's{'m'^)]c is the contribution of the branch cut. 
The contribution of the pole is given by 



7^ dx dpi.p\ exp(^-PY^p2 _^ „^2^ 



(2^) 

Fi{m'^,-m'^)y/V + F2{m'^,-m'^)Qs{i\lpi+rn?^)^ 
X ^ — ^ (A.29) 

+ m?^/V 
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with l'^ = U jV'^l'^ and -p^ = p±V^/^. The contribution of the branch cut is given by 



[JUm% 



dx I d^p± / dpo exp{iL^po) 

J i^yp^+m^ /x{l—x) 

DiscJ^l(m^p^) + Disc J^2(?n^ff^)gs(poV'^/^) 
p"^ + 



(A.30) 



where 



Disc J^i(m2,/) = lim [Ti{pI + (po + ef) - Ti{p\ + (po - e)^)] 

£-5-0 '- -' 

= 2-Ki[w? + x{l — x)p'^\ , 
Disc J2(m2,p2) = [j-2(p^ + (po + e)2) _ ^2(^1 + {po - ef)] 

£-5-0 '- -' 

= Airix'^ . 



(A.31) 



Therefore 



(2^)2y 



I dx / #_L]5i / 

JO J,/)3' 



dyexpi-l^y) 



m 



^p\+m?\/V /x{l-x) 

V + - x){p\ - y^) + 2x^gsiiy)VV 

y2 _ p2^_ ffl^^^Y 



(A.32) 



with Po = iy and thus dpo = idy. In Sec. 3 we calculate [Pg^*" (m2)]finite numerically at scale 
/i = i.e., we replace Fi{vn? ^p^) and F2{m?' ^p^) by 



(A.33) 



J-i{m^ ,p^) = [m^ + x(l — log[m^\/y + x{l — x)p^\/v'\ 
T2{m^,P^) = 2x^{l + \og[rr?^/V + x(l - x)p^^/V]) . 

A.2.3 The term Pgi^(m2) 

The term Pl^{m?) is equal to the term P^^ {iv?). This can be seen by shifting the integration 
variables p^ Pfi — and using the invariance of the integral under A;^ — )• — fc^. 

A.2.4 The term P^^{m'^) 

The term Pg^^(m2) is finite and can be calculated numerically. We rewrite P^^{m'^) analogously 

to P^^^im^) as 



dxdydz 



{xy + xz + yzy 



X exp 



X exp 



2{x + y) 

E 



{Ll/VV)[-2soy + 2roix + y)]^7T 



yz + xy + xz 
{L]/^/V)[-2rjSjy + r|(x + y) + s]{y + z)]?!' 



{x + y + z) 



47r 



(A.34) 
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In the following we discuss how to efficiently calculate the infinite sums over sq, • • • , S3, tq, . . . , r^. 
We first define 



9= Yl n ^M-ajr] - ujoj 



bjs'j +CjrjSj +djrj + ejSj) 



= n exp(-ajr| - bjs] + cjrjSj + djrj + ejSj) - J| exp(-6js| + ejSj) 

s^O ^ r j 

- n exp[-(aj + bj - Cj)s] + (dj + ej)sj] 
j 

= ^ + n exp(-ajr| - 6jS^ + cjrjSj + djTj + ejSj) -^Yl ^M-bjs'j + e^sj) 

r,s j s j 

- ^ JJexp(-ajr| + djrj) - ^ exp[-(aj + 6^ - Cj)s'j + (dj + ej)sj] , (A.35) 
i « j 

where aj , 6j , cj , dj , Cj £ C with Re aj , Re 5j > and j = 0, . . . , 3. In this way we can write 

3 3 3 3 

j=0 i=0 j=Q j=0 

(A.36) 

where 

z^(a, b) = exp(— an^ + bn) , 

n 

5(0, 6, c, d, e) = exp(— an^ — bm? + cmn + dn + em) 

n,m 

= ^{a, cm + d) exp(— + em) . 

m 

Since is a Jacobi theta function it transforms covariantly under inversion of a, 



(A.37) 



z^(a, 6) = W — exp(6 /4a)i^(7r /a,iii{b/a)), 



which is readily shown using Poisson's summation formula. For a G IR, with a < vr we can 
therefore use Eq. (A.38) to achieve a swift convergence of the sum over n in v{a, b) in Eq. (A.37). 
In fact in the least favorable case of a = vr we need only to sum all n with |n — ?io| < 4, where 
n = no G Z maximizes —av? + hn, to achieve a precision of — log|o[exp(— 4^7r)] k, 22 digits. In 
this way we can express Pg^^(m^) as 



Pl'^im') 



(Svr)- JO 
X 5(60,^0) 



dxdydz R 



eo=do=0 



exp 



2(x + y)- Rli[-2yde, + 2{x + y)dd,Y7r 



(A.39) 





(i) 


(ii) 


(iii) 


extrapolation 


(oi) 


f\ r\-t rw cvi 

—0.010591 


—0.010538 


r\ r\'\ r\f f\'\ 

—0.010501 


—0.0106(1) 


(«3/2) 


0.009395 


0.009429 


0.009439 


0.00940(4) 


(02) 


0.024867 


0.024892 


0.024884 


0.02487(2) 




U.Uoiy44 


n HQ 1 von 


n HQ 1 I'VQ 

U.Uoi / 10 


U.Uoiy^i j 


(04) 


-0.021362 


-0.021531 


-0.021720 


-0.0214(3) 


(^3/2) 


-0.015129 


-0.015092 


-0.015050 


-0.01513(3) 


(^2) 


-0.015765 


-0.015714 


-0.015669 


-0.01577(7) 


(^3) 


-0.008005 


-0.007968 


-0.007904 


-0.0080(1) 


(^4) 


0.011724 


0.011748 


0.011810 


0.0117(1) 



Table 3: Best fit of Pg for regions (i), (ii), and (iii) and extrapolated values of Pg for different geometries. 



with 



R = , Ij = L/v^/^ , 

yz + xy + xz 



Rlj {x + y)TT , bj = Rlj {y + z)-k , 



Rl']2-Ky , di = d2 = ds = ei = 62 = 63 = . (A.40) 



The integral over x, y and z can be performed conveniently in spherical coordinates, as in the case 
A.3 The complete diagram 

We combine all contributions to PQ{m'^) and find that the complete diagram at scale fi = V'^l^ is 
given by 

(-P6)uv =3-^ — 3~'^'^^ ' 

- + 2[P6''^(m2)]fi„,e - + Pl''\m') + Pi^im') , 

P^ = lim P^{m^) , (A.41) 

where Pq is the renormalized shape coefficient at scale n = V~^l'^. In Figs. 4 and 5 we show the 
result of a numerical calculation of Pg (yn^) for different values of rr?\fV . 

Note that the linear- divergences as well as the logarithmic divergences in 1/ rr? VV cancel. 
We perform a fit to a polynomial of order four for different ranges ofm^^/V■. (i) m'^^/V < 0.75, 
(ii) m^W < 1.5, and (iii) rn^VV < 3. The result of the fits and the corresponding extrapolated 
values for = are given in Table 3. 
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Figure 4: Extrapolation of Pg = lim„j^o ("^^) ^'-'^ lattice geometries {ax), x G {1, 3/2, 2, 3, 4} at scale 
/i = We fit a polynomial of order four including numerical data from three different ranges (i), (ii), 

and (iii). 
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ih) (64) 



Figure 5: Extrapolation of Pg = lim„i_j.o -Pe (™^) for lattice geometries (fe^), x E {3/2, 2, 3, 4} at scale 
fi = V^^^'^. We fit a polynomial of order four including numerical data from three different ranges (i), (ii), 
and (iii). 
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